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ABSTRACT 

We report the discovery of accretion disks associated with ~ 13 Myr-old 
intermediate/low-mass stars in h and x Persei. Optical spectroscopy of ~ 5000 
stars in these clusters and a surrounding halo population reveal 32 A-K stars with 
H Q emission. Matching these stars with 2MASS and optical photometry yields 
25 stars with the highest probability of cluster membership and EW(H Q ) > 5 A. 
Sixteen of these sources have EW(H Q )> lOA. The population of accreting sources 
is strongly spectral type dependent: H Q emission characteristic of accretion, es- 
pecially strong accretion (EW(H a ) > 10A), is much more prevalent around stars 
later than GO. Strong H Q emission from accretion is typically associated with 
redder K s -[8] colors. The existence of accreting pre-main sequence stars in h and 
X Persei implies that circumstellar gas in some systems, especially those with 
primaries later than G5 spectral type, can last longer than 10-15 Myr. 

Subject headings: Galaxy: Open Clusters and associations: Individual: NGC 
Number: NGC 869, Galaxy: Open Clusters and associations: Individual: NGC 
Number: NGC 884, accretion, accretion disks, planetary systems: protoplanetary 
disks 

1. Introduction 

Young stars are born with massive, ~ 0.01-0.1M* (stellar mass), disks of gas and dust. 
The disk viscously spreads, transporting angular momentum away from the star and mass 
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onto the star. Accretion onto the host star is identified from strong H a emis sion; typical mass 
accretion rates onto the star are ~ 10~ 8 M yr _1 for ~ 1 Myr-old stars (IHartmann et al 



19981 ). After ~ 5 Myr, fewer sources show strong H a emission indicative of accretion, and 
accretion rates are typically much lower (~ 10~ 9 yr -1 ) than at earlier ages. By ~ 10 
Myr, few sources show signs of active accretion ( Sicilia-Aguilar et al. IboOol ). 



The timescale for accretion to cease and for nebular gas in circumstellar disks to dis- 
perse has important i mplications for planet formation. Massive planets may f orm in gas- 
poor/free conditions (ICurrie and Hansen 1 120071 ; iKonacki and Wolszczan 1 120031 ). However, 
gas giant planet formation requires that circumstellar gas lasts longer than the time for 
a ~ 10-15 Mp protoplanetary core to capture a large gaseous envelope (~ 10 6 -10 7 years; 



Ikoma et al. Il2000l ; iPapaloizou and Nelson 



shorter, the disk instability model (e.g. iBoss 1 120051 ) for Jovian planet formation also re- 



2005). Though the formation timescale is much 



quires a massive gas disk. Circumstellar ga s is also necessary for planetary migration (see 
Goldreich and Tremaine Ill980t IWard 1119971 ). which may explain close-in giant planets ('hot 
Jupiters'), and circularization of the orbits of terrestrial planets. A spread in the time for 
accretion to cease and for gas to disperse may then lead to diverse planetary systems. 



While previous studies show that few stars accrete gas after ~ 10 Myr fjSicilia-Aguilar et al 



20051 ). the bulk properties of the longest-lived accreting systems are not well constrained 
due to small number statistics. There are ~ 5 known actively accreting stars with ages 
> 10 Myr, including PDS 66 AK Sco, and St 34 (Mamajek et al. 2002; Reipurth et al. 
1996; Iwhite and Hillenbrand Ihoosf ). Well- known 10-20 Myr old clusters such as NGC 7160 
Jsicilia-Aguilar et al. Ilioiil ) and Sco-Cen (IChen et al. 1120051) have < 200 know n members. 
Because accretion at ages > 10 Myr is rare (~ 2%; ISicilia-Aguilar et al. 1 120051 ) . a larger (> 
1000-2000 sources) ensemble of 10-20 Myr sources spanning a wide range of spectral types 
is required to investigate the oldest accretion disks. 



The double cluster, h and y Persei - 13 ± 1 Myr old, d=2.34 kpc (jBragg and Kenyon 



20051 ; ISlesnick et al.ll2002l ; Keller et al.ll200ll ) - provides an excellent opportunity to study 
the properties of a statistically significant sample of the longest -lived accretion disk s. With 
> 5000 members in the cluster and surrounding halo population (ICurrie et al.ll2007al ). robust 
constraints on gas accretion are possible even if the fraction of accreting sources is ~ 1-2%. 

In this Letter, we report the discovery of 25 sources in h and y Persei with evidence for 
active gas accretion. In §2 we identify these accreting sources by measuring H a equivalent 
widths and use quantitative spectral types to select actively accreting h and y Per members 
from a sample of ~ 5000 optical spectra. Analysis of our sample in §3 reveals that the 
population of accretors is almost solely comprised of stars with spectral types later than 
G5. We compare the strength of accretion signatures to dust emission in §4, finding a 
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weak correlation between the H Q equivalent width and IR excess. We conclude with a brief 
summary and discuss future observations to investigate circumstellar gas accretion in young 
stars in more detail. These results show that circumstellar gas in at least stars later than ~ 
GO lasts > 10-15 Myr. 



Hectospec Observations/Hydra Archival Data of h and \ Persei and 

Sample Identification 



We obtained Hectospec (IFabricant et al.l 120051 ) spectra of 4536 sources with V ~ 16-19, 
J~ 14.25-16.25, and J-H ~ 0-1.5 in x Persei and the halo population surrounding both clus- 
ters on the 6.5m MMT telescope at Whipple Observatory during September and November 
2006. Each source was observed in three 10-minute exposures using the 270 mm -1 grating. 
This configuration yields spectra at 4000-9000 A with 3 A resolution. The d ata were pro- 
cessed using the standard Hectospec reduction pipeline (IFabricant et al.ll2005l ) and typically 
had S/N > 30-50 at 5000 A. 

We acqui red additional spectra of 710 sources near h and \ Per with the Hydra multifiber 



spectrograph (IBarden et al. 



19931 ) on the WIYN 3.5 m telescope at the Kitt Peak National 



Observatory. Hydra spectra were obtained during two observing runs in November 2000 and 
October 2001 and include stars with V ~ 14-17 and J ~ 12-14.5 (S/N « 10-30 typically). 
We used the 400 g mm -1 setting blazed at 42°, with a resolution of 7 A and a coverage of 
3600-6700 A. The standard IRAF task dohydra was used to reduce the spectra. 



To select candidate emission line stars, we measured spectral indices (1(A); see lBalog and Kenyon 
2002h with the IRAF sbandE routine. Because classical T Tauri stars have strong H a emis- 



H 7 , H,5, and several 



sion and often have strong Ca II emission, we derived indices for H a , H 
Ca II features. For the H I features, bandpasses with widths of 30 A ( H a ) and 20 A (Hp, 
H 7 , and H^) yield good results. For Ca II, we centered bandpasses at 3933 A (20 A width) 
and at 8498, 8542, and 8662 A (20 A) for the IR triplet lines. To identify sources with H a 
emission, we compared the spectral index of H a to Hg, which was less likely to be in emission 
and typically had high signal-to- noise. 

The H Q and Hg spectral indices for most stars lie in a narrow band on I(Hg) vs. I(H a ) 
from ~ (-0.2,0) to (0.7,0.4). Acc reting sources have anom olo usly small H n index for a 
given Hp index (Figure dh; see also lBalog and Kenyon 1 120021 and lBragg and Kenyon 1120021 ). 



1 IRAF is distributed by the National Optical Astronomy Observatory which is operated by the Associ- 
ation of Universities for Research in Astronomy Inc. under contract to the National Science Foundation. 
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Candidate accreting sources are those lying outside the main distribution of sources on I(Hg) 
vs. I(H a ), below the line I(H Q )=0.5x I(H J g)-0.15 (Figured^). After removing background 
M giants, Be stars, and low signal-to-noise stars, 32 sources with clear H Q emission and 
moderate to high signal-to-noise spectra remain. 

While ea ch source clearly shows Hq, emission, sources with weaker emission may not 
be accreting (IWhite and Basri 1 120031 ) . Chromospheric activity can produce very small H a 
equivalent widths (EW (H a )< 3-5 A). Accreting T Tauri stars typically EW(H a ) > 10-20 
A, though low rates of accretion (M > 10~ 8 M^yr" 1 ) can produce EW(H„) ~ 3-10 A for 



stars earlier than K5 (IWhite and Basri 1120031 ; iBarrado y Navascues and Martin 1 120051 ) . We 
consider sources with EW(H a ) > 5A as potential accreting stars. To distinguish more clearly 
between accreting sources and chromospherically active sources, we remeasured EW(H a ) 
using the IRAF routine splot, smoothing any low signal-to-noise spectra and fitting the line 
profile to a gaussian. Typically, splot derived only slighly smaller EW(H a ) than sbands (~ 
a few A). Sources with EW(H a ) > 10A are 'strong' accret ors. Consistent with observat ions 
of T Tauri stars in Taurus-Auriga and other clouds (e.g. Kenyon and HartmannlE995l ). no 
sources show evidence of Call emission. 

To derive quantitative spectral types for these 32 sources, we used spectral indices for 
Hp, if 7 , Hs, the G band (4305 A), and Mg I (5175 A). W e fit piecewise l i near r elationships 



between each spectral index and spectral type from the iJacoby et al. I (119841 ) standards, 
derived spectral types for each index, and adopted the median as the spectral type. For 
the 3 sources with Hp in emission, we derived a spectral type from the 4 other indices 
without signs of accretion. The resulting spectral type distribution ranged from A7 to M0 
with typical uncertainties of ~ 1-2 sub classes. We compa r ed ou r results to those using the 



quantitative spectral typing method of iHernandez et al. 
to be excellent. 



(120041 ) and found the agreement 



To identify stars likely associated with th e clusters, we used o ptical photometry (IKeller et al 
200lh and 2MASS/IRAC photometry from ICurrie et al.l (b()07ah . We exploited the well- 



constrained age and reddening of the double cluster and derived the expected 2MASS J 
magnitude of members at each spectral ty pe. Sources with in 0.2 mags of a band defined by 
the 12 and 14 Myr stellar isochrones (from lsiess et alihoool ) were identified as members (see 
Currie et al.ll2007al ). Figure [T] shows the 30 sources consistent with membership: the other 
two sources were a late K star and a M0 star that are too cool to fall within the isochrone. Al- 
though these two sources may lie at the cluster distance and have slightly diffe rent ages (e.g. 



10 or 15 Myr) or extinctions, the clusters appear to have a small age spread (iMeynet et al. 



19931 : IKeller et al.ll200ll ; ISlesnick et al]l200a Bragg and Kenvon 1120051 ) . Thus, the 30 sources 
shown in Figure [lb are the most likely cluster members. Table 1 lists their properties. 
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Y oung T Tauri stars show a st rong correlation between EW(H a ) and IR excess from 
disks (IKenyon and Hartmannl Il995l ). Stars in ~ 1 Myr-old Taurus-Auriga with L band 
(3.8/im) excess (K-L > 0.4) typically have stronger H Q emission (EW(Ha) > 10A) than 
those lacking excess. Excess sources have a range of EW(H Q ), which implies a range of 
accretion rates. To search for a similar trend in o ur sample, we use the K s -[8] color derived 
from the 2MASS/IRAC data (ICurrie et al.ll2007al ) as an indicator of excess emission instead 
of K s -[3.6]. Eighteen of 30 sources with H Q emission are detected at [8] with errors < 0.2 
mags. Typically, sources without 8/im photometry are fainter in K s (~ 14.25-14.75) than 
those with S^im detections (K s ~ 13.5-14). 



3. Results: Evidence for Accretion at ~ 13 Myr 

Most of the H Q emission stars in h and x P er have GO or later spectral types (Figure 
H]). Overall, 25/30 sources have EW>5 A and thus are likely accreting circumstellar gas. 
All fourteen sources later than G5 have EW >5 A. By contrast, 5/16 sources G5 and earlier 
have EW(H a ) < 5A . The distribution of sources with EW(H Q ) > 10A is even more strongly 
biased towards later spectral types. Seven out of eight sources later than K0 have EW(H Q ) > 
10A while only 6/16 of the sources G5 and earlier have EW(H a ) > lOA. Four of five sources 
with EW(H Q )> 20A are later than K0, and the sources with the strongest H Q emission (EW 
~ 47 A and 59 A) are a K7 and K4 star. The atlas of both emission line and non-emission line 
stars will appear in the full spectroscopy survey paper (T. Currie et al. 2008, in progress). 

The spectral type distribution of accreting sources is comprised almost exclusively of 
sources later than GO. Because we do not have a complete sample of non-accreting sources at 
a given spectral type, we cannot accurately measure the frequency of accretors as a function 
of spectral type/mass. Nevertheless, we can estimate the fraction of accreting sources in our 
sample by deriving how many of the 5145 (4536 from Hectospec, 609 from H ydra) sources 



with ne ar-IR photometry lie along the h and x P er isochrone in J/J-H. We follow lCurrie et al. 



(l2007al ) and identify as cluster members 2040 sources within a band that extends 0.3 mags 
fainter and 0.75 mags brighter (the binary locus) than the 13 Myr isochrone, where the vast 
majority of non members lie below the isochrone. We estimate the spectral types of sources 
from their J-H colors and divide the sample into two populations, earlier than F8 and later 
than F8, and compute the frequency of accreting sources in each population. 

The frequency of accreting sources is low and is consistent with a spectral-type depen- 
dence. Overall, ~ 1.2% (25/2040) of sources consistent with cluster membership are likely 
accreting. Through F8, 1/1042 (0.1%) sources with spectra show evidence of accretion. The 
frequency is 2.4% (24/998) for later sources. Spectral types for the non-accreting sources 
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will provide better constraints on the frequency of long-lived accretion disks. The low frac- 
tion of accreting stars in h and \ Per is broadly consistent w ith results from other studies. 
Sicilia-Aguilar et al. I (120051 ) and ISicilia-Aguilar et al. I (120061 ) identified accretion in ~ 40% 
of stars in the 4 Myr-old Tr 37 Cepheus OB2 subgroup and only one of 55 stars (~ 2%) in 
the other Cep OB2 subgroup, NGC 7160 (~ 11.8 Myr old). An accretion disk fraction of ~ 
1.2% from our sample is consistent with a decline of accretion disk frequency with age. 

To check these results, we consider possible selection biases. If the emission-line stars 
have significant optical veiling, their spectral types appear earlier than fainter (J < 16) 
non-accreting stars with similar masses. We then overestimate their relative frequency. 



. „ , „ — ^„ 

(Kenvon and Hartmann 


1995) 


2007a 


), emission-line stars in 



negligible veiling. If photometric errors move faint cluster stars outside the isochrone, we 
underestimate the number of non-accreting stars. Including all stars with J-H > 0.48 (F8; 
2303 total) in the cluster sample yields a strict lower limit to the accretion frequency, 1.1% ± 
0.2%, which is still much larger than the frequency for earlier type stars. Finally, we consider 
the possibility that the Ha emission flux associated with late-type stars is undetectable 
around 1 mag brighter cluster stars with earlier spectral types. With EW(Ha) ~ 5-10 A 
for the late-type stars, early-type stars with similar Ha fluxes should have EW(Ha) = 2-5 
A, easily detectable on our spectra. Thus, we conclude that possible selection effects do not 
produce a false trend of increasing accretion frequency among lower mass stars. 

Our sample shows a ~ 2-3<x correlation between accretion and infrared excess. The 
majority of stars have negligible K, 9 -[8] excess (F i gure [3]) , but some sources have K,,-[8l~ 0.4 , 
the nominal cutoff for IR excess in ICurrie et all (I2007ah and iKenyon and Hartmann! (119951 ) . 
One source, the K7 star, has K s -[8] ~ 1 and is clearly redder than a stellar photosphere. 
Sources with larger H a (>20A) have redder K s -[8] colors than those with small/marginal 
H a emission. The Spearman rank correlation coefficient between K s -[8] and H Q emission is 
r s =0.60; this distribution has a low probability of being drawn from a random sample (p,i 
= 0.9%). The distribution is less correlated for sources with chromospheric H Q emission 
(r s =0.35, pa ~ 22%). We thus find a ~ 2-3 a correlation between IR e xcess and EW(H fV ) f or 
sources with H a emission. However, some debris disks in h and x P er ( Currie et al. 2007bl lcl) 
have strong excesses at [8] and/or [24] without accretion. Therefore, while accretion may 
imply IR excess, IR excess need not imply accretion. 
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4. Summary & Discussion 



We have discovered a population of stars with accretion disks in h and x Persei. Of the 
30 H a emitters that are likely h and x Per cluster/halo members, 25 have EW(H Q ) consistent 
with circumstellar gas accretion. The spectral types of accreting sources are almost all later 
than GO. About ~ 1.2% of our sample show s evidence for accretio n , cons istent with the low 
fraction of accretors at ~ 10 Myr found by ISicilia-Aguilar et al. I (120051 ) . We find that H a 
emission sources with larger EW(H a ) have redder K s -[8] excess. 

While current survey points to a spectral type dependence on the accretion disk lifetime, 
a full survey of h and x Persei sources can provide stronger constraints on the properties 
of long-lived disks. Our sample is drawn heavily from x Persei and the surrounding halo 
population. Analyzing spectra of the slightly more massive/dense member of the Double 
Cluster, h Persei, will test whether the frequency of accreting stars is the same or depends 
on properties such as the mass and density of the cluster. A complete Hectospec survey 
of ~ 15,000 sources in the h and x Per field is underway and will better characterize the 
population of accreting pre-main sequence stars. 



We thank the anonymous referee for a timely and helpful report and Eric Mamajek 
for useful discussions regarding previously known stars with long-lived accretion disks. TC 
and SK acknowledge support from the NASA Astrophysics Theory Program grant NAG5- 
13278 and TPF Grant NNG06GH25G. ZB received support from Hungarian OTKA Grants 
TS049872 and T049082. 



REFERENCES 

Barden, S., et al., 1993, PASP, 37, Fiber Optics in Astronomy II, ed. P. M. Gray (San 
Francisco: ASP), 185 

Barrado y Navascues, D., Martin, E., 2005, RMxAC, 24, 179 

Balog, Z., Kenyon, S., 2002, AJ, 124, 2083 

Boss, A., 2005, ApJ, 629, 535 

Bragg, A., Kenyon, S., 2002, AJ, 124, 3289 

Bragg, A. & Kenyon, S., 2005, AJ, 130, 134 

Chen, C, et al., 2005, ApJ, 623, 493 



- 8- 

Currie, T., Hansen, B., 2007, ApJ, 666,1232 

Currie, T., et al., 2007(a),ApJ, 659, 599 

Currie, T., et al., 2007(b), ApJL, 663, 105 

Currie, T., et al, 2007(c), ApJ accepted. la7Xiv:0709.2510l 

Fabricant, D., et al., 2005, PASP, 117, 1411 

Goldreich, P., Tremaine, S., 1980, ApJ, 241, 425 

Hartmann, L., et al., 1998, ApJ, 495, 385 

Hernandez, J., et al, 2004,AJ, 127, 1682 

Ikoma, M., et al., 2000,ApJ, 537, 1013 

Jacoby, G., et al, 1984, ApJS, 56, 257 

Keller, S.C., et al, 2001,AJ, 122, 248 

Kenyon, S., Hartmann, L., 1995, ApJS, 101, 117 

Konacki, M., Wolszczan, A., 2003, ApJ, 589, 495 

Mamajek, E., Meyer, M., Liebert, J., 2002, AJ, 124, 1670 

Meynet, G., Mermilliod, J-C, Maeder,A., 1993, A&AS, 98, 477 

Papaloizou, J., Nelson, R., 2005, A&A, 433, 247 

Reipurth, B., Pedrosa, A., and Lago, M., 1996, A&AS, 120, 229 

Siess, L., et al, 2000,A&A, 358, 593 

Sicilia-Aguilar, A., et al., 2005, A J 

Sicilia-Aguilar, A., et al., 2006, A J 

Slesnick, C, et al., 2002,ApJ, 576, 880 

Ward, W., 1997, Icarus, 126, 261 

White, R., Basri, G., 2003, ApJ, 582, 1109 

White, R., Hillenbrand, L., 2005, ApJ, 621, 65 



- 9- 



This preprint was prepared with the AAS IAT^X macros v5.2. 



-10- 



Table 1. H a emission sources 



a 8 ST cr(ST, subclasses) EqW(H Q ) Accreting? J K s K g -[3.6] K„-[8] 



2:22:21.62 57:04:00.5 G4 1.0 


13.7 


y 


14.30313.612 0.292 




2:22:39.53 57:15:42.9 GO 2.5 


4.5 


n 


14.31513.538 0.271 


0.223 


2:17:36.96 57:01:10.6 GO 3.0 


1.8 


n 


14.36313.697... 


-0.247 


2:17:55.15 56:56:03.1 G3 2.1 


6.46 


y 


14.41913.613 0.258 


0.096 


2:21:46.44 57:02:44.2 G4 2.3 


5.75 


y 


14.44613.6970.349 


0.405 



Note. — Properties of sources with H a emission (first five entries). Sources with EW(H a ), 
> 5A are identified as accreting sources. All but the final entry are from Hectospec. Un- 
certainties in EW(H Q ) were ~ 1-2 AAll photometric measurements listed are 5 a detections. 
Sources were typically brighter than the 10c limit in J (15.7) and K s (14.8) (see lCurrie et all 
l2007ah . 




Fig. 1. — (left) The H Q spectral index plotted against the Hg index. Sources lying below the 
dotted line were selected for emission signatures; sources between the main distribution and 
the line typically had lower signal-to-noise, (right) Spectral type vs. J magnitude for h and 
X Per members with H Q emission. The dotted line indicates the faint limit of the survey. 
The 12 and 14 Myr isochrones (top and bottom thick lines) are overplotted with 0.2 mag 
errors (top and bottom grey lines). The errors in derived spectral type are overplotted as 
bars on the data points (triangles). 
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Fig. 2. — Spectral type distribution of H a emission sources (grey line), those with EW>5A 
(dash-dot line), and those with EW>10A (thick dotted line). Sources with H a ~5-10A are 
probably marginally accreting, those with H a > 10A are consistent with the strength of H a 
emission from Classical T Tauri stars. 
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Fig. 3. — H Q equivalent width vs. K s -[8] color. There is a weak correlation between the H, 
width and IR excess. 



